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1 INTRODUCTION 



Isospin violation in decays of e^e~ B^B^ at the T{AS) resonance results in a difference between 

the branching fractions /oo = B{e^e~ B^B^) and /_| = B{e^e~ — > B~^B~). Measurements of 

the ratio R^^'^ = /_| //ooj summarized in Tabled are consistent with unity within the errors pQ. 

Theoretical predictions for i?"*"/*^ range from 1.03 to 1.25 [2]. Currently, almost all published mea- 
surements of B meson branching fractions make the assumption that R^/^ = 1. Precision mea- 
surements of /oo, /h , and R^/^ can be used to eliminate this assumption and re-normalize all B 

meson branching fractions. 



Table 1: Summary of previous measurements of R'^/^ 



Decay B 


JCdt 




Source 


J/HK^/Kl!) 


81.9 fb-i 


1.006 ±0.036 ±0.031 


BABM 3 




2.73 fb^i 


1.058 ±0.084 ±0.136 


CLEO g 


J/^/i(+/o) 


20.7 fb-i 


1.10 ±0.06 ±0.05 


BABAR 




9.2 fb-i 


1.04 ±0.07 ±0.04 


CLEO i 



In this paper we report the first direct measurement of /oo- The measurement is based on partial 
reconstruction of the decay B^ D*'^i~Ui (the inclusion of charge-conjugate states is implied 
throughout this paper). This allows a sizeable sample of double tagged events to be identified. 
Comparison of the double-tag and the single-tag yields allows a determination of /qq with minimal 
input from simulation. 

The technique used to measure /oo is as follows: in every event we reconstruct the decay B^ 
D*~^£~v'i, as described further below. The sample of events in which at least one B^ — > D*~^i~i'i 
candidate decay is found is labeled as "single-tag sample" . The number of signal decays found in 
this sample is 

Ns = 2iV^^/oo es B{W ^ D*+ri?e), (1) 

where = (88726 ± 23) x 10^ is the total number of BB events in the data sample and is 
the reconstruction efficiency of the decay B^ D*~^i^i'£. The technique for measuring is 
described in [7\. The data sample has a mean energy of 10.580 GeV IE,; and an energy spread of 
only 4.6 MeV . Such a small spread means that any energy dependence of /oo has a negligible effect 
on the central value. The subset of single-tag events in which two B^ D*^i~Ui candidates are 
found is labeled as "double-tag sample" . The number of such events is 

Nd = /oo [BiB'' ^ D*+rD,)]\ (2) 

where is the efficiency to reconstruct two B^ —>■ D*^£~i'£ decays in the same event. Note that 
every double-tag event contributes two entries to the single-tag sample. Using Eq. (Q) and Eq. 
the ratio /oo is given by 

where we have defined the coefficient C = ed/e^- C = 1 if the efficiencies for detecting each B 
meson are uncorrelated in double-tag events. 
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2 THE BABAR DETECTOR AND DATASET 



The BABAR data sample used in this paper consists of 81.7 fb^ cohected at the T(4S') resonance 
(the on-resonance sample) and 9.6 fb~^ collected 40 MeV below the resonance (the off-resonance 
sample). Simulated BB events were analyzed through the same analysis chain as the data. The 
equivalent luminosity of the simulated sample is approximately three times that of the on-resonance 
data. 

A detailed description of the BABAR detector and the algorithms used for track reconstruc- 
tion and particle identification is provided elsewhere [H]. A brief summary is given here. High- 
momentum particles are reconstructed by matching hits in the silicon vertex tracker (SVT) with 
track elements in the drift chamber (DCH). Lower momentum tracks, which do not leave signals 
on many wires in the DCH due to the bending induced by a magnetic field, are reconstructed by 
the SVT alone. Electrons are identified with the ratio of the track momentum to the associated 
energy deposited in the calorimeter (EMC), the transverse profile of the shower, the energy loss in 
the drift chamber, and the information from a Cherenkov detector (DIRC). Muons are identified in 
the instrumented flux return (IFR), composed of resistive plate chambers and layers of iron. Muon 
candidates are required to have a path length and hit distribution in the instrumented flux return 
and energy deposition in the EMC consistent with that expected for a minimum-ionizing particle. 
The Cherenkov light emission in the DIRC is then employed to further reject kaons misidentified 
as muons by requiring muon candidates to have a kaon hypothesis probability less than 5%. 

Hadronic events are selected by requiring at least four charged particle tracks reconstructed 
by the silicon vertex detector and the drift chamber. To reduce background from continuum 
e~^e~ — > qq, where q stands for a u, d, s, or c quark, the ratio R2 = H2/HQ of the second to the 
zeroth Fox- Wolfram moments is used |inj . 

3 ANALYSIS METHOD 

We reconstruct the decays B^ D*^£~U£ with a partial reconstruction technique. The application 
of the technique to this mode was first proposed by the ARGUS Collaboration and has been 
used by CLEO DELPHI _12j, OPAL and BABAR il4|. In this technique, only the lepton 
from the decay B^ D*^£~i>^ and the soft pion from the decay D*^ —>■ D^tt~^ are used. No 
attempt is made to reconstruct the D^, resulting in high reconstruction efficiency. 

To suppress leptons from charm decays, all lepton candidates (electrons and muons) are required 
to have momentum between 1.5 GeV/c and 2.5 GeV/c in the e~^e~ center-of-mass (CM) frame. 
Soft pion candidates are required to have CM momentum between 60 MeV/c and 200 MeV/ c. As 
a consequence of the limited phase space available in the D*~^ decay, the soft pion is emitted 
within a one radian-wide cone centered about the D*^ direction in the CM frame. The D*~^ four- 
momentum can therefore be computed by approximating its direction as that of the soft pion, and 
parameterizing its momentum as a linear function of the soft pion momentum, with parameters 
obtained from the simulation. The presence of an undetected neutrino is inferred from conservation 
of momentum and energy. The neutrino invariant mass squared is calculated: 

= (^beam " ^D* - E,f - (pz,* + , (4) 

where E-^^^^ is the beam energy and {Ed*) and (pd*) are the CM energy and momentum 
of the lepton (the D* meson). If the decay is properly reconstructed and the neutrino is the only 
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missing particle, the distribution will peak near zero for signal events. Background events, 
however, are spread over a wide range of Ai'^ values. 

In what follows, we use the symbol to denote M'^ for any candidate in the single-tag sample. 
In the double-tag sample, we randomly choose one of the two reconstructed 0*^1" vi candi- 

dates as "first" and the other as "second". Their M"^ values are labeled Ad\ and A^l; respectively. 
For each of the variables ^A1 {i = s,l,2), we define a signal region M.f > —2 GeV^/c^ and the 
sideband -8 < < -4 GeVVcl 

In addition to signal — > D*'^£~ui decays, the single-tag and double-tag samples contain 
several types of events: 

• Continuum e+e~ qq background. 

• Combinatorial BB background, formed from random combinations of reconstructed leptons 
and soft pions. This background can also be due to the low momentum soft pions not coming 
from a D*, produced by either the same B or other B [T5] . 

• Peaking BB background, composed of i? ^ D* (mr)ir'£ decays with or without an excited 
charmed resonance (D**) [TBI, where the reconstructed soft pion comes from the decay D*~^ — >■ 
D^7r~^, leading to an accumulation of these events at high values of Mf. These events are 
peaking background and are produced both by B^ and B~^ decays. Their Mf distribution 
differs from the signal, which allows us to extract their contribution in a fit. Such events are 
suppressed by the requirement pi > 1.5GeV/c on the lepton CM momentum. 

• The decays B^ — > D*^T~i'r and B^ D*^ i~ r'i[n'y) may be used for the measurement of /oo 
and are therefore considered as signal. These events peak in A4f and come from B^ decays. 

The double-tag sample contains two additional types of background: A^f -combinatorial and Al^- 
peaking. In Alf -combinatorial (A4f-peaking) background events the first candidate is combinatorial 
(peaking) background. 

To determine Ns and N^, we perform binned fits to one-dimensional histograms of the A4'^ 
and distributions of on-resonance data events, ranging from -8 to 2 GeVVc^- Before fitting, 
we subtract the continuum background contribution from the histograms. This is done using the 
A4l and Al^ distributions of off-resonance data, scaled to account for the ratio of on-resonance 
to off-resonance luminosities and the CM energy dependence of the continuum production cross- 
section. In addition, the contributions of the Alf -combinatorial and A^f-peaking backgrounds are 
subtracted from the A^| histogram before the fit. The A1 ^-combinatorial background is determined 
from the A4f sideband, which contain only continuum and combinatorial background events. This 
histogram is scaled by the ratio of the number of combinatorial events in the signal region and the 
sideband, determined from the simulation. The Al^-peaking background subtraction is based on 
the simulated Alf-peaking events. 

After the subtraction, the A4^ and AI2 histograms are fit separately, using a function whose 
value for bin j of the histogram is 

f^=Y^N'pj, (5) 

t 

where A^* is the number of events of type t {t = signal, combinatorial, peaking) populating the 
histogram, and Pj is the bin j value of a discrete probability density function (PDF) obtained from 
simulated events of type t, normalized such that Pj = 1. The fit determines the parameters A^* 
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by minimizing: 

y2 = V ~ ^^^^ (6) 

3 tj 



where Hj is the number of entries in bin j of the histogram being fit; is the statistical error 
on Hj, including uncertainties due to the background subtractions described above; and af. is the 
/j, determined from the errors on Pj. 



error on /,-, determined from the errors on which are due to the finite size of the simulated 



sample. 

The results of the fits are presented in Table El The and M.2 distributions are shown in 
Fig. n with the contributions of the different event types indicated. The fits yield the values Ng = 
786300 ± 1950 and Nd = 3560 ±80. Using the simulation we determine C = 0.9946 ± 0.0078, where 
the error is due to the finite size of the simulated sample. Eq. ^ then gives /oo = 0.486 it 0.010, 
where the error is due to data statistics only. 



Table 2: Numbers of entries of different types found by the fits to the M.1 and histograms 
in the signal region. Also shown are the numbers of entries of subtracted backgrounds and the 
confidence levels of the fits. 



Source 


Mi 


Mi 


Combinatorial BB 


558090 ± 760 


1520 ± 40 


Peaking BB 


68170 ± 260 


300 ± 20 


Signal 


786300 ± 2000 


3560 ± 80 


Continuum 


238500 ± 1300 


160 ± 40 


A^f -combinatorial 




180 ± 20 


Alf-peaking 




60± 10 


xVd.o.f. 


41/56 


48/56 


Confidence level 


93% 


77% 



To determine how well the simulation reproduces the M.1 and M.^ distributions of the combina- 
torial background in the data, we study the distributions of a sample of same-charge candidates, in 
which the lepton and soft pion have the same electric charge. This sample contains only continuum 
and combinatorial BB background. We fit the continuum-subtracted M.1 and M.^ histograms of 
the same-charge sample using the function /j = NP'^, where Pj is the bin j value of the PDF of 
same-charge simulated BB events, normalized such that ~ ^1 parameter is deter- 

mined by the fit. The histograms, overlaid with the fit function, are shown in Fig. |5J The ratio 
between these two histograms is fitted to a constant both for the M.^ and M^ summed over the 
signal region and over all bins are shown in Fig. |3J The accumulated differences D = ^j{Hj — /j) 
between the same-charge data histograms Hj and the fit functions are summarized in Table 01 Their 
consistency with zero indicates that the distributions of simulated combinatorial BB background 
events do not lead to significant fake signal yields. 

4 SYSTEMATIC STUDIES 

We consider several sources of systematic uncertainties in /qq. All estimated errors are an absolute 
systematic uncertainties in /oo and summarized in Tabled 



11 



Table 3: The difference D = J2jif^j ~ fj) between the same-charge data histogram and the fit 
function, summed over all bins or over the signal region only. Also shown are the fit x^/d-o.f. 
values and confidence levels. 





Signal rej 


;ion 


All bins 






Ml 


-Mi 


Ml 


Ml 


D 


-1300 ±2100 


-80 ± 80 


-700 ± 3000 - 


-70 ± 80 


xVd.o.f. 


17/19 


13/19 


40/55 


34/53 


C.L.(%) 


59 


84 


94 


98 



Table 4: Summary of the absolute systematic errors for /qq. 



Source 




A^^-combinatorial 


0.0005 


A^f -peaking 


0.0005 


Same charged events 


0.0025 


Peaking background 


0.004 


i?-meson counting 


0.0055 


r(45) uou-BB 


0.0025 


Efficiency correlation 


0.004 


Monte Carlo statistics 


0.002 


Total 


0.009 



1. The systematic uncertainty from the 7W|-combinatorial contribution subtraction in the Ml 
histogram is 0.0005. The error is obtained by varying the total A^f -combinatorial background 
by its statistical error and repeating the analysis. 

2. An error of 0.0005 is estimated due to the subtraction of the Alf-peaking contribution in 
the Ml histogram. The error is obtained by comparing the ratio between the numbers 
of subtracted A4f-peaking and A4f -combinatorial events with their ratio of peaking and 
combinatorial events in Table |2l 

3. Propagating the errors on the quantities D (same charged events) of Table El leads to an error 
of 0.0025 on /oo- To determine this error we vary the signal events both for the single-tag 
and the double-tag samples. The largest uncertainty then is taken for the uncertainty on /oo- 

4. The PDFs P* (t = peak) of the peaking background come from simulated event samples 
containing different D** resonances or non-resonant events. We vary the ratio of the branching 
fraction of the resonant and the non-resonant production such that the variation of this ratio is 
wide enough to include poorly known decays. We repeat the analysis procedure to determine 
Ng and A''^. The resulting error on /oo is 0.004. 

5. Uncertainties in the branching fractions of D*^t~Vt- and B^ D*~^ l~ u^{n'^) relative 
to B^ D*^i~r'£ lead to uncertainties in the PDFs P* {t = signal) of the signal events. This 
uncertainty in /oo is negligible. 
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6. The error due to the uncertainty in is 0.0055. It includes the uncertainties for differ- 
ences in the cross sections and efficiencies for muon pairs and continuum events between 
on-resonance and off-resonance samples, hadronic selection criteria and the uncertainties in 
the tracking efficiency. 

7. In this paper the impact of non-BB decays of the T(45) on i?-meson counting has been 
accounted for as a systematic error. The upper limit for the branching fraction of T{4S) 
decays into non-BB is 4% at 95% confidence level We conservatively estimate the 
systematic error by decreasing the 4% upper limit on the branching fraction to 2%. From 
this variation we estimate an error of 0.0025 due to the effect on of a possible decay. 

8. We note that the lepton momentum spectrum in the Monte Carlo simulation is different from 
the one we observe in the data. We tune the simulation to the data by rejecting simulated 
events in such a way that the two lepton momentum spectra agree. We repeat the analysis 
procedure without the rejected events. The systematic error due to the uncertainty in the 
lepton momentum spectrum is negligible. 

9. There is a small efficiency correlation between the single-tag and the double-tag samples. 
The systematic uncertainty due to this efficiency correlation is estimated by propagating the 
Monte Carlo simulation systematics error of C into /oo- The simulation statistical error in C 
leads to a 0.004 error in /qq. In addition to the Monte Carlo simulation systematics error of 
C, we study the effect of track multiplicity on the efficiency correlation. 

10. We perform a similar procedure as mentioned above for the pion momentum spectrum. The 
error due to the uncertainty in the pion momentum spectrum is negligible. 

11. An error of 0.002 is due to the finite size of the simulated sample, calculated using af . in 
Eq. ©. 

12. The estimator used in Eq. ® can be biased. We did an alternative binned likelihood fit 
and found that the result differed by only 0.03% for /qq. 

We combine the uncertainties given above in quadrature to determine an absolute systematic 
error of 0.009 in /qq. 

5 SUMMARY 

To summarize, using partial reconstruction of the decay B^ D*^l^i>i we have obtained a pre- 
liminary result for the branching fraction 

/oo = 0.486 ± 0.010 ± 0.009, (7) 

where the first error is statistical and the second is systematic. Since this measurement is done by 
comparing the numbers of events with one and two reconstructed B^ — > D*~^(.~vi decays, it does not 
depend on branching fractions of the B^ and the D*^ decay chains, on the simulated reconstruction 
efficiency, on the ratio of the charged and neutral B meson lifetimes, nor on assumptions of isospin 
symmetry. 
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Figure 1: The Ai^ (top) and (bottom) distributions of the on-resonance samples. The con- 
tinuum background has been subtracted from the distribution. For the M2 distribution, the 

A^f-combinatorial, and the TWf-pcaking have been subtracted. The levels of the simulated signal, 
peaking BB and combinatorial BB background contributions are obtained from the fit. 
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Figure 2: The Ai'^ (top) and Al^ (bottom) distributions of the same-charge on-resonance samples. 
The continuum background has been subtracted from the distribution. For the A^l distri- 
bution, the continuum background, the A^f-combinatorial and the A^f-peaking backgrounds have 
been subtracted. The level of the simulated combinatorial BB background is obtained from the fit. 
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Figure 3: The ratio between data and the combinatorial BB background of the same-charge sample 
both for the Ai'^ and summed over the signal region and over all bins. The values are fit to a 
constant, upper left: for the Ad^ summed over the signal region; upper right: for the summed 
over the signal region; lower left: for the A4g summed over all bins; lower right: for the summed 
over all bins. 
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